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We study the interplay between an inhomogeneous quantum quench of the external potential in 
a system of relativistic fermions in one dimension and the well-known Klein tunneling. We find 
that the large time evolution is characterized by particle production at a constant rate which we 
derive analytically. The produced particles can be physically interpreted according to a semiclassical 
picture and the state reached in the long time limit can be classified as a non-equilibrium-steady- 
state. Such a quantum quench can be used in order to observe macroscopic effects of Klein tunneling 
in transport, with potential implementations with current experimental setups. 
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I. INTRODUCTION 

Klein tunneling p] refers to the fact that an incoming 
relativistic electron can penetrate a potential barrier of 
height V greater than twice the electron’s rest mass m if 
the electron’s energy lies between m and V — m 1 . In par¬ 
ticular, keeping the energy of the incident particle fixed 
and increasing the height of the barrier, the transmis¬ 
sion coefficient does not decay, as one would probably 
expect and as it actually happens in the bosonic ana¬ 
logue [2j. The key feature of this phenomenon is the fact 
that the electrons are relativistic and in fact the neces¬ 
sary tool to correctly explain it is quantum field theory 
m\- After its explanation (see mmm for more recent 
reviews), the Klein tunneling phenomenon gained much 
relevance along with the discovery of particle-antiparticle 
production from a strong enough potential [5], vacuum 
polarization effects [9] and black hole evaporation due to 
the creation of particle-antiparticle pairs near the event 
horizon (Hawking radiation) [lOj . 

This and other peculiar effects of the Dirac equation, 
like the Zitterbewegung EL although representing key 
phenomena to understand relativistic quantum effects, 
have proven difficult to observe experimentally. For in¬ 
stance, the observation of Klein tunneling requires a po¬ 
tential drop of the order of the fermion mass m over the 
Compton length 1/m which yields an enormous electric 
field [T2] '13], thus making the effect relevant only for 
very exotic situations mm- These difficulties stimu¬ 
lated a great interest for the simulations of relativistic 
quantum systems with condensed matter setups in the 
laboratory |TS| . The pioneering study of graphene [H HU 
is certainly the primary example. More in general, us¬ 
ing recently developed techniques, like ultracold atoms 
in optical lattices [H] , ions E] or photonic systems [20], 


1 Throughout the paper we will use natural units h = c = 1. 


it is possible to device highly tunable systems which al¬ 
low for the preparation and detection of a great variety 
of many-body phenomena. Examples include the sim¬ 
ulation of black holes in Bose-Einstein condensates [21] , 
Dirac equation in various dimensions |22| . Recently Klein 
tunneling was also simulated in single ion traps 123] . 

In this paper we study another physical situation where 
the features of Klein tunneling emerge via a sudden 
quench of the external potential for a system of relativis¬ 
tic one-dimensional Dirac fermions from a constant zero 
value to a step-like profile. 

The sudden quench protocol [23] consists of the fol¬ 
lowing three steps: (1) prepare a system in a pure state, 
usually the ground state of a pre-quench Hamiltonian; 
(2) at a certain time, namely t = 0, suddenly vary one 
of the parameters of the Hamiltonian; (3) from that time 
on, let the system evolve unitarily, i.e. without connec¬ 
tion or dissipation to the external environment, accord¬ 
ing to the post-quench Hamiltonian which does not com¬ 
mute with the pre-quench one. Being the pre-quench 
and post-quench Hamiltonian not mutually commuting 
guarantees that the time evolution starts from an initial 
out-of-equilibrium configuration. 

This is the simplest way of driving a many body quan¬ 
tum system out-of-equilibrium and has the advantage 
that analytical calculations of the long time limit of ob¬ 
servables can be carried out in many cases [25] Ej. When 
the pre-quench Hamiltonian is inhomogeneous, it is also 
possible to study transport properties of the system [27]. 
It is important to highlight that many of these recent 
theoretical investigations are triggered by several exper¬ 
imental realizations of quantum quenches [28] in engi¬ 
neered low dimensional quantum systems in ultracold 
atoms [T3]> which allow to have coherent dynamics for 
much longer times than with usual solid state systems. 

The reason why Klein tunneling can be linked to the 
quantum quench scenario relies on the characterization 
of the initial state in which the system is prepared before 
the quench in terms of the post-quench Hamiltonian. It 
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is a well-established semiclassical picture [2U (2Hj that by 
performing a global quench on the system (i.e. changing 
some parameter in the whole system) we inject an exten¬ 
sive amount of energy and the initial state can be thought 
of as containing an extensive number of quasi-particle ex¬ 
citations which, after the quench, start spreading ballisti- 
cally. In the present case we show that, as a consequence 
of quenching the potential in the one dimensional Dirac 
equation from a constant zero value to a step-like profile 
with height V > 2m, at large times there is particle pro¬ 
duction at a constant rate, which we derive analytically. 
Our rigorous field theoretical derivation also provides a 
solid theoretical framework for previous numerical inves¬ 
tigations of particle production [Ml El] and single particle 
quantum mechanical calculations JB] that perfectly agree 
with our findings. 

The manuscript is organized as follows: in section ( |TT| ) 
and (III) we carefully re-derive the energy eigenstates 
that are solutions of the Dirac equation in a homoge¬ 
neous potential and in an inhomogeneous one of a step¬ 
like form, pointing out the differences depending on the 
height of the potential. In section (IV) we introduce the 
quench of the potential from V = 0 to V(x) = V0(a:), 
in particular we derive the overlaps between the pre¬ 
quench and post-quench eigenfunctions; then we analyt¬ 
ically compute the number of pre-quench particles pro¬ 
duced through the scattering process at large times and 
give a physical interpretation to it. 


II. ONE-DIMENSIONAL DIRAC FERMIONS IN 
PRESENCE OF POTENTIAL STEP 

A. Homogeneous case 

We consider a one-dimensional system of relativis¬ 
tic fermions of mass m in an external potential V ( x ). 
This is described by the Dirac equation ( i<J) — to — 
7 o V(x))t/j(x,t) = 0 , which, using the representation 
7 o = cr 3 , 71 = zcri, reads 

((id t - V(x))a 3 - (7 1 d x - m)ip(x, t) = 0, (1) 

where < 71 , 2,3 are the Pauli matrices. The energy eigen¬ 
states satisfy the time-independent form of the Dirac 
equation 

[(E - V(x))a 3 - aid x - m]i/)(E-,x) = 0, (2) 

where ip(E-,x) is an eigenstate of energy E. 

In the absence of a potential, i.e. when V(x) = 0, the 
solutions to the above equation are 

= «*(£;*) (3) 

where fc® = VE 2 — to 2 with |i5| > m. Positive energy 
solutions E > +m describe particles while negative en¬ 
ergy solutions E < —to describe antiparticles. For each 


energy E > +m (E < — to) there are two independent 
solutions corresponding to left or right moving particles 
(right or left moving antiparticles, respectively). Note 
that the above solutions are orthonormalized in infinite 
volume i.e. 

+ OO 

J dx {u s (E-x)) ] u s ' {E'-x) =S ss> S(E-E') (4) 

— OO 

If the potential V (x) is homogeneous i.e. V(x) = Vo, then 
the solutions are ip{E\x) = ^(E — Vq;x)-, consequently 
particles are described by solutions with E > V +m while 
antiparticles by solutions with E < V — m (see Fig. [I]). 


particles 


m 


—m 

anti-particles 


V 


V + m 


particles 


V — m 


anti-particles 


FIG. 1: Energy windows for the eigenfunctions in the ho¬ 
mogeneous case, with distinction between the particle and 
antiparticle solutions. On the left, the ones for the V = 0 
case, on the right the ones for the V case. 


B. Inhomogeneous case 

Assuming now a potential step V(x) = VQ(x) (with 
V > 0) the new solutions can be found by matching to¬ 
gether at the origin the homogeneous solutions in the 
positive and the negative semi-axes with the correspond¬ 
ing values of the potential, keeping in mind that solutions 
with imaginary wavenumbers are also acceptable in the 
semi-axes as long as they decay exponentially at large 
distances. The matching is prescribed by the continuity 
condition ip(E\ 0~) = i/j(E;0 + ) at the origin. Depending 
on the value of E , the wavenumbers q = ks-v on Vhe 
left and p = kE on the right may be both real, one real 
and the other imaginary or both imaginary. We there¬ 
fore distinguish the following cases, which are depicted 
for clarity in Fig. [2] 

1) V < 2m: In this case there exist three different en¬ 
ergy windows. For energies E £ (— 00 , —m] U[V +m, + 00 ) 
both q and p are real and the solutions correspond to 
plane waves on both sides of the step. They repre¬ 
sent particles or antiparticles that cross the potential 
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V <2 m 


—m V — in m V + m 



V >2 m 


-m m V -m V + m 

m 



FIG. 2: Energy windows for the different kinds of eigenfunc¬ 
tions in the inhomogeneous potential case. Grey zones denote 
totally reflecting solutions, green ones plane wave solutions. 
The crossed one, only present in the V < 2m case, represents 
the not allowed region, while the dashed one, present only in 
the V > 2m case, represents the Klein Zone. 


(a) 



iP) 



E £ (—m, V — to) U (to, V + m) instead, exactly one of 
q and p is real while the other is imaginary, therefore 
the solutions correspond to plane waves on one side of 
the step and exponentially damped waves on the other. 
They represent particles or antiparticles that are totally 
reflected by the potential step and are non-degenerate 
since, in order to have exponential damping on one side, 
the incident plane wave can come only from the other 
side. Lastly, for energies Eg [V — m, to] both q and p 
are imaginary and so the associated solutions, if there 
existed, would be bound states. However the continuity 
condition is not satisfied for any value of the energy in 
this window, so there do not exist any solutions of this 
type. 

2) V > 2m: In this case the energy windows are as 
follows: solutions are plane waves on both sides for ener¬ 
gies E £ (—oo, —to] U [to, V — to] U [V + to, +oo), while 
for all other energies E £ (—to, to) U (V — to, V + m) the 
solutions are totally reflecting. 

Note that, as already known from standard quantum 
field theory, in the ground state the infinite set of negative 
energy eigenstates should be considered as occupied up 
to some energy level, the Fermi sea level Ep, in order for 
the energy spectrum to be bounded from below, so that 
the theory makes physical sense. Excitations above the 
ground state are either occupied eigenstates with energy 
E > Ep or unoccupied eigenstates with energy E < Ep 
(hole or antiparticle excitations). Ep is typically chosen 
to be zero, however in the present problem we would 
rather let ourselves free for the moment. 

Let us first focus on the case V < 2m. The (doubly 
degenerate) reflecting-transmitting solutions are given by 


vi(E;x) 


v 2 {E\x) 


u + (E;x) + Aiu (E;x), x<0 
B 1 u + (E — V;x), x>0 


(5) 


B 2 U (E; x), x < 0 

u~(E — V; x) + A 2 U + (E — V; x), x>0 

( 6 ) 


where, using the matching condition, the coefficients A ± t 2 
and B 12 are found to be 


. 1 - k 

A 1 ~ 1 , ! 

1 + K 

A2 = —Ai, 


B 1 

B 2 


2y/~K 
1 + K ’ 

Bi, 


(7) 

( 8 ) 


with 


q {E + to) 
p (E — V + to) ’ 


(9) 


step and are partially reflected and partially transmit¬ 
ted to the opposite side. Since the incident particle 
or antiparticle may come from either side of the ori¬ 
gin, the solutions are doubly degenerate. For energies 


and 

q = k E -v = V(E - V ) 2 - to 2 , (10) 

p = kE = \/ E 2 — to 2 . (11) 
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The reflection R and transmission T probability coef¬ 
ficients for both V\ .2 are given by 


R = \A\ 2 = 


(12) 

T = \B\ 2 = 

4 M 

(1 + k) 2 ' 

(13) 


For energies E G (— 00 , — m] U [V + to, + 00 ), the param¬ 
eter k is positive, hence Ai and Bi are real and therefore 
R + T = 1. This relation expresses the conservation of 
probability currents. 

The physical meaning of ( [5| a nd © for E > V + m is ex¬ 
plained respectively in Fig.j3]a) and (b). For E > V + m 
the eigenfunctions represent particles both in (|5j) and ©> 
so the group velocity v g = dE d ^ , indicated by the arrow, 
is in the same direction as the momentum (indicated by 
the index ± of the u function). In Fig. [3ja) , for x < 0 
u + and u~ are respectively the incoming (from the left) 
and reflected (to the left) particles, while for x > 0 u + is 
the transmitted (to the right) particle. In Fig. [3](b) , for 
x < 0, u~ is the transmitted (to the left) particle, while 
for x > 0 u~ and u + are respectively the incoming (from 
the right) and reflected (to the right) particle. 

The totally reflecting solutions are given by the same 
expressions, more precisely by ([ 5 ]) if E G (to, V + m) and 
by © if E G (—to, V — to) with the coefficients A t and 13, 
still given by (18]) but now being complex. Ai in particu¬ 
lar is unitary (i.e. |A;| 2 = 1) since n is imaginary. The 
reflection and transmission coefficients are in this case 
R = 1 and T = 0. Note that these are non-degenerate 
solutions: the first solution corresponds to a particle in¬ 
cident from the left and exists for E G (to, V + to), 
while the second corresponds to an antiparticle inci¬ 
dent from the right (with coefficient A 2 ) and exists for 
E G (—to, V — to). The non-degeneracy is because 
it is only one of the two complex-conjugate imaginary 
wavenumbers ±k,E-v = ±i\ltE-v\ (or +ks = ±z||) 
that results in a wavefunction that is exponentially de¬ 
caying at large distances x —> +00 (or x —> —00 respec¬ 
tively) . 

It can be verified that the above eigenstates are by con¬ 
struction orthonormalized, since vA 1 are orthonormalized 
too. In particular the pair of degenerate states ^ and 
([6]) is chosen in such a way that they are orthogonal to 
each other 

+ OO 

J dx (vx(E-,x))^ v 2 (E'-x) = 

— OO 

\ {A\B 2 + BIA 2 )5(E - E') = 0, (14) 

since A\ and B\ are real. Normalization is ensured by 
the probability conservation relation R + T = 1, 

+OO 

J dx (v s (E;x))^ v s (E';x) = 

— OO 

\{l + R + T)6{E-E')=6(E-E'), s = 1,2. (15) 


III. KLEIN TUNNELING 


We now turn our attention to the case V > 2m. While 
it is still true that for energies E G (— 00 , —to] U [V + 
to, + 00 ), the parameter k is positive and R+T — 1, this is 
no longer true within the energy window E G (to, V — to), 
which we will call from now on ‘Klein zone This is 
because the parameter k becomes negative and so R and 
T as defined in (12), (13), now satisfy R — T =1. The 


probability conservation seems then to be violated since 
R + T 7 ^ 1. Moreover increasing the potential step V 
the parameter k, for energies in the middle of the zone 
(to, V — to), tends to -1 so that T tends to unity. 

As first clarified by Sauter [3] and Hund [J, this phe¬ 
nomenon is paradoxical only as long as it is understood in 
a one-particle quantum mechanical framework but there 
is no room for any paradox within second quantized 
field theory, after taking into account the possibility of 
particle-antiparticle production When an external 

potential is high enough it allows for the spontaneous pro¬ 
duction of particles-antiparticles pairs □HE]; the bar¬ 
rier is repulsive for electrons but attractive for positrons. 
Then Klein tunneling arises when fermion states outside 
the barrier have energy that matches the energy of the 
antiparticles inside the barrier imm- Since electrons in 
x < 0 have energies larger than to while positrons in 
x > 0 have energies smaller than V — to, a non van¬ 
ishing overlap between these two regions is present only 
for V > 2to and actually represents the Klein zone (see 
Fig. [4]. 

As a consequence, for all energy zones where k > 0 the 
eigenfunctions maintain the same form as in the case 
V < 2m. On the other hand, for V > 2m and k < 0, 
the coefficients Ai are real but Bi are imaginary, so that 
neither orthogonality (14) nor (15) holds for the states 
©and§. ‘ • LJ LJ 

The correct ones can be found noticing that, for the 
state V\ to describe a scattering process with a transmit¬ 
ted (i.e. right-moving) antiparticle on the right, the di¬ 
rection of momentum on that side must be flipped. Sim¬ 
ilarly, for v 2 both directions of the plane waves on the 
right should be flipped (in which way the amplitude of 
the incoming particle is still equal to 1). We thus define, 
in the Klein zone, the new pair of degenerate states 

= ( 16 ) 

L B[u (E—V;x), x > 0 


v 2 (E-,x) = 


B' 2 u (E'x), x<0 

u + (E — V; x) + A' 2 u~ (E — V] x ), x > 0 

(17) 


and using the matching condition we find that the new 
coefficients A[ 2 and B[ 2 are 


, 1 + n 1 — k' 

1 1 — k 1 + k' ! 

A' — - A’ 


= = (18) 


1 — K 1 + k’ 


B’ 2 = -B[, 


(19) 
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a.)V <2m b) V > 2m 




FIG. 4: Absence (a) and presence (b) of Klein tunneling respectively for the cases V < 2m and V > 2m. We highlight the 
energy window for particles in x < 0 (light-blue), for antiparticles in x > 0 (green). Only for V > 2m there is a non vanishing 
overlap between the two (the blue region). 


where 


q (E + m) 
p (E — V + to) 


> 0 . 


( 20 ) 


Now the reflection R and transmission T probability 
coefficients, defined as the reflected-to-incoming and 
transmitted-to-incoming probability current ratios, are 
given by 


R=\A'\ 2 = 

\1 + k'J ’ 

(21) 

T = \B'\ 2 = 

I 

(1 + K') 2 ’ 

(22) 


for energies in the Klein zone, 


and since k' > 0, the probability conservation relation 
R + T = 1 is recovered. 

The new states (16) and ( [T7| ) are orthogonal to each 
other 


+oo 

J dx (v[(E;x)^ v' 2 (E'-x) = 


±(A'* 1 B' 2 + B'* 1 A' 2 )6(E-E') = 0, 


(23) 


after taking into account (19) and that B[ is imaginary. 
Each of these states is normalized to unit 


+oo 

J dx {v' s {E-,x)) ] v' s (E'-,x) = 

— OO 

§(1 + R + T)5{E-E') = 5(E-E'), s = l,2. 

(24) 


It should be stressed that redefining the eigenstates 
for the Klein zone as above is a necessity rather than 
an arbitrary choice: indeed, for the subsequent study or¬ 
thogonality and normalization are necessary. If we chose 
to keep one of the eigenstates, say r; 2 , in the same form 
as in the other zones, we would have to multiply it by 
1 / A 2 in order to make it normalized in the Klein zone, 
thus obtaining v 2 = V2/A2, and next to choose the other 
eigenstate as the (unique) linear combination of v\ and 
V 2 that is orthogonal to v 2 and normalized too, thus ob¬ 
taining = v\ — (i?i / A 2 )v2 ■ The new pair of states v[ 
and v 2 are then precisely the ones given by (161 and (JTT 


Notice that since the group velocity of an antiparticle is 
opposite to its momentum, all eigenstates corresponding 
to antiparticles should be defined, for any value of the 
potential, with their momentum signs flipped, if we want 
them to describe physical scattering processes in which 
the incident wave has coefficient equal to 1. When all 
three plane waves (incoming, reflected and transmitted) 
correspond to antiparticles (for instance in the case E < 
—m of ([5j) and Q) flipping all three momentum signs of 
the eigenstates does not amend the values and physical 
significance of their coefficients A*, Bi or those of R and 
T, neither does it spoil the orthonormalization relation 
of the two degenerate states. On the contrary, in the 
Klein zone, since in x < 0 there are particles while in 
x > 0 antiparticles, only the momenta of the plane waves 
in x > 0 must be flipped. For this reason, while it is 
crucial to redefine the states in the Klein Zone, this is 
not necessary for all other energy windows and we can 
keep them in the original form © and Hence, as 
a pictorial representation for the eigenfunctions ([16 1 and 
© in the Klein zone, we can refer to Fig. [ 3 ] with the 
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index of u in x > 0 changed. The arrows represent the 
correct direction of the group velocity of plane waves. 

Except when otherwise stated, in the following we will 
drop the prime from the notation of expressions (161 and 
that is we redefine 'ty to be equal to v[ in the Klein 
zone energy range and similarly for the parameters Ai 
and B;. 


IV. QUENCH 

We now consider a one-dimensional system of Dirac 
fermions described by initially prepared in the 
ground state that corresponds to a homogeneous poten¬ 
tial Vq = 0; at t = 0 the potential is quenched to the 
inhomogeneous step potential V(x) = V@(x) and the 
system, from now isolated from the rest, is subject to 
unitary evolution according to the post-quench Hamilto¬ 
nian. 

The pre-quenclr Hamiltonian in a second quantized 
form is 

H = E ( r dEEa E a E- [ m dEEp$/3%), 

a =± \Jm J— oo / 

(25) 

where a E and f3 E are the electron annihilation opera¬ 
tor and positron creation operator respectively, obeying 
usual anticommutation relation 


{a E , a:^, ] 

■ = $a,<j'8{E — E') t 

{ a E , a E'} ~ 0, 

(26) 


>=5 a ^5(E-E'), 

'Co 

II 

0 

(27) 

The post-quench Hamiltonian is 



H = E ( 

[ dE E <4 s)t E - 

- [ dEEb^b^ 

). 

EE 

J £ > 

' £ < 

J 

(28) 

where 




We, a E'} 

\=5 s ^S(E-E'), 

{ a Ei a E'} = 0; 

(29) 

{b s E ,b%} 

[ = 5 S ^5{E-E'), 

ft* 

O- 

II 

O 

(30) 


An excitation that corresponds to a particle occupying 
an eigenstate v s (E;x) with energy E > E F = 0 is asso- 
ciated with a creation operator a E , while an excitation 
that corresponds to an unoccupied eigenstate v s (E;x) 
with energy E < Ep = Oisa hole excitation and is 
associated with a hole creation operator b y E . In the 
above, £> and £< denote the energy ranges over which 
there exist solutions of type s = 1,2, above or below 
the considered ground state level Ep respectively. More 
explicitly, as can be seen in Fig. [2j eigenstates v\ ex¬ 
ist for energies E £ (—oo, — to) U (+m, +oo), while i >2 


exist for E £ (—oo, V — to ) U (V + m, + 00 ). For ener¬ 
gies E £ (V — m, V + to) the state v\ becomes totally 
reflecting but they are still given by the same form. Sim¬ 
ilarly, for energies E £ (—to, +m) the state V 2 becomes 
totally reflecting. Lastly, for energies E £ (m, V — m) the 
states vi } 2 are given by (16) and (p~7|). Ep is chosen to 


be the one corresponding to the pre-quench Hamiltonian 
i.e. Ep = Vo = 0. Overall we have £> = [+m,+ 00 ), 
£> = [0, V — to] U [V -f- m,+ 00 ), £< = (— 00 ,— to] 
and £3 = (— 00 ,0]. The Klein zone will be denoted 
by £ KZ = [m,V — m]. Similarly, we will denote by 
£> = (+m, + 00 ) and £< = (— 00 , —to) the energy ranges 
over which there exist pre-quench solutions u a (a = ±) 
above and below Ep = 0. 


We are interested in the evolution of the total number 
of pre-quench particles, i.e. electrons, at large times after 
the quench, when considering the thermodynamic limit 
of our system ( i.e. the limit L —> 00 with L the length 
of the system). We would like to point out that, being 
the initial state of the whole system a pure one and the 
time evolution unitary, the system will always be in a 
pure state at any time, exhibiting quantum recurrences 
in its time evolution which, in the present problem, are 
due to particles moving around the circumference L of 
the system. On the other hand, first taking the thermo¬ 
dynamic limit and then the long time limit, finite sub¬ 
systems of the whole system can be described by a mixed 
state , whose observables usually exhibit a stationary 
limit [33, 34]. This prescription is the standard way of 
deriving the large time asymptotics in extended quantum 
systems. It is then clear that the only relevant observable 
to our analysis is the pre-quench number of particles, not 
the post-quench number of excitations, whose time evo¬ 
lution is trivial. 

The evolution of the fermionic field i!){x,t) is formally 
given by expanding on the post-quench creation and an¬ 
nihilation operators. 


ip{x,t) = 



+ j dE v s (E-x)b { E }t 
■£< 


e~ iEt 


(31) 


In order to calculate expectation values of observables 
after the quench, we need to know the initial expectation 
values of the post-quench creation and annihilation oper¬ 
ators and these can be found expressing them in terms of 
pre-quench creation and annihilation operators, by com¬ 
parison of the expansion of the field ip(x,0) in the two 
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different bases. In the pre-quench basis we have 


^0,0) = 


E 

O'—± 


+oo — m 

J dE u a (E; x)a E + j dE u a (E\ x)/3^ 


E [ dEu G {E-x) 7 |. 

a=± J 


(32) 


where we defined the operator 



\ a E 

if E > m, 


tqs 

III 


if E < —to, 

(33) 



otherwise. 


Clearly, a E and /3^ are the electron annihilation 

opera- 


tor and positron creation operator respectively. 
Using the orthonolmalization relations, we find that 


+oo 

J dx (u a {E-,x)]* ij){x, 0 )=Je- (34) 


Similarly, in the post-quench basis we have 
ip(x, 0) = 


E 

s—1,2 


dE v s {E\x)a' E + / dE v s (E;x)b^ 


s= 1,2 

where 


/ dE v s (E;x)c 

,— i o 


.(») 

E ■ 


4 s) 

4 S) = <1 & ( E s)t 


if Ee 
if £*, 
otherwise. 


(35) 


(36) 


The orthonolmalization relations lead now to 


are the overlaps between pre-quench and post-quench 
eigenstates. For convenience and brevity, we will incor¬ 
porate the energy zone limits into the expressions for the 
overlaps, defining 

w as {E, E') = W as (E, &)$£. (E)d S s (E’) (40) 


where the function ds{E) equals unit if E £ £ and zero 
otherwise. 

Similarly the inverse of (38) is 



(41) 


and expresses the post-quench operators in terms of the 
pre-quench ones. 


A. Overlaps between pre-quench and post-quench 
eigenstates 


We will now explore the properties of the overlaps 
W as (E,E') that will be essential for the subsequent cal¬ 
culation. Substituting the expressions for the eigenstates 
and taking into account that 


+oo 


dx e lkx O(x) 


lim -- 

e->o+ k + ie 


(42) 


we find that the overlap functions have poles whenever 
the two energies match so as kE = fcg' ± ie or k E = 
—kE> ±ie. Explicitly, W±\{E, E') for E' outside the Klein 
zone is given by 


W±i(E,E') 


iD ± (E, E') 

— ks 1 + *e 


iD^(E,E')A 1 {E') + iD ± (E, E' - V)B 1 {E') 
ikE + kE' + + kE'—v + ie 

(43) 


+oo 

J dx (v s (E-x))^ i/j{x, 0) = Ce\ (37) 


From (34) and (351 we can write the pre-quench annihi¬ 


lation and creation operators in terms of the post-quench 
ones 


1e= E / dE'W„ s (E,E')c$ 

s=l 2/, 


(38) 


where 


+ OO 

W as {E,E') = J dxu^{E-x)v s (E'-,x) (39) 


and for E' in the Klein zone 

iD ± (E,E') 


W ±1 (E,E') = 


EkE — kE’ + ie 
i J CU(F;,£’ , )Al(F; , ) W T (E, E' - V)B 1 {E’) 


±k E + ks ' + ie 


— kE’-v + ie 


(44) 


where 

D ± {E,E’) = 


1 

4-7T 


E + to E' + m 


kE 


kE 1 


± 


kE 


kE' 


E + m V E' + m 


(45) 


























Similarly the overlaps W± 2 {E, E') for E' outside the 
Klein zone are 


W ±2 (E, E') 


iD T (E, E' - V) 
Tks — kE'-v + 


iD ± (E, E' - V)A 2 {E') 
+ kE<-v + 


iD*{E,E')B 2 (E') 

±fc.E + kE> + 

(46) 


and for E' in the Klein zone 


w± 2 (E, E') 


iD ± {E,E' - V ) 
-fkE + kE'-V + 


iD*{E,E' - V)A 2 {E') 
=F/c_e — kE'-v + 


iD*{E,E')B 2 (E') 

± ks + kE 1 + 

(47) 


Written as functions of the energy E, the overlaps 
have simple poles close to the real axis at E ^ E' and 
E « E' — V . All other poles and branch-cut singular¬ 
ities that are located away from the real E -axis do not 
matter in the thermodynamic limit, since their contribu¬ 
tion is exponentially suppressed in the L —> oo, as can 
be seen from the fact that e ~ 1 /L (see Appendix [AJ . In 
order to express the overlaps around the above poles, it 
is sufficient to substitute 


ip~\E') 


kE^ks’+ie <Je'{E — E') + iep 1 (E I ) 
ip-\E') 


+ 


-a E '{E + E')+iep^{E') 


for real kE,kE' 


(48) 


where <je is the sign of E , p(E) = \E/ks | is the den¬ 
sity of states at the energy E and the dots denote 
corrections that are functions non-singular along the real 
E-axis. Terms proportional to i/(kE + kE> + ie), which 
do not have poles close to this axis can be omitted. More¬ 
over, it turns out that there are no poles at opposite ener¬ 
gies E ss —E' or E « —(E' — V), because the correspond¬ 
ing residues are proportional to E + (E, —E) = 0. This 
expresses the fact that particles and antiparticles with 
the same absolute value of energy are orthogonal. The 
residues of the poles at E ss E' or E ~ E' — V are pro¬ 
portional to E + (E,E) = p(E)/(2n) times the amplitude 
of the incoming, reflected or transmitted wave. Terms 
proportional to i/(kE ± fcs' + ie) with E' £ (—m,+m) 
which correspond to energies in the totally reflecting 
zones, can be omitted since kE 1 is imaginary in this en¬ 
ergy window, so it cannot match with which is al¬ 
ways real. 

After some algebra, and taking the above substitution 
rules into account, we find that in all but the Klein zone, 
the overlaps can be written in the form 


W +a (E,E') 


W_ i(E,E') 


._ voo _, 

a- e{E - E') + ie/p(E) 
iB 1 (E+V)/{ 2tt) 

-a E {E -E' + V) + ie/p{E) + 

iA 1 (E)/(2n) 

-<j e {E - E') + ie/p(E) 
iA 2 (E + V)/(2it) 

- (j e {E -E'+V) + ie/p(E) + 


(49) 

(50) 

(51) 


W_, 2 (E,E') 


_V(2tt)_ 

' a e (E-E' + V)+ie/p{E) 
tB 2 (E)/(2n) 
-a E (E-E')+ie/p(E) 
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(52) 


while in the Klein zone 


W +:1 (E,E’) 


W_,i {E,E') 


W +t2 (E,E') 


W_, 2 (E,E') 


._ i/V*) _ , 

<t E (E - E’) + ie/p(E) 
iA[(E)/(2n) 

-<j e {E - E') + ie/p{E) 
iB[(E + V)/(27t) 

(t e {E -E' + V) + ie/p(E) 

._ VOO _ , 

■ -<je(E -E' + V) + ie/p(E) 
iA! 2 (E + V)/( 2 tt) 

‘ (t e {E - E' + V)+ ie/p(E) 
iB’ 2 (E)/( 2n) 

-cte(E — E') + ie/p(E) 


(53) 


(54) 

(55) 


(56) 


where, as above, “...” denote corrections that do not in¬ 
volve any poles on or close to the real .E-axis and primed 
quantities have been used to emphasise that we refer to 
the Klein zone expressions. Note that the residues Ai,Bi 
and A ! i , B[ are smooth and bounded functions (except 
at the edges of the energy zones, where they are still 
bounded, but discontinuous or non-smooth). The above 
poles express the resonance that occurs when the in¬ 
coming, reflected or transmitted wave of the post-quench 
eigenstate has the same energy and direction as the pre¬ 
quench eigenstate wave and is of the same type (i.e. par¬ 
ticle or antiparticle). 


B. Evolution of the number of particles 


We now proceed to the calculation of the evo¬ 
lution of physical observables. We will focus on 
the total number of pre-quench particles N(t) = 
T,a=±J^ dE and, in particular, its 

behaviour in the large time and thermodynamic limit. 


W +t2 (E,E') = 
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OO 

m = E /dE m)hE^E\m) = 

°=±L 

OO 

E [ dE / ^1^2 < ai (£?,£!)«;„,(£?,£7 2 ) ($(0)|4 1 i t 4 2 J$(0)) = 

<7,s i ,S2 m 

OO 

E l dE J dE.dEzdE' w* asi (E,E 1 )w aS 2 (E,E 2 )w <Tlsi (E / ,E 1 )wl, S 2 (E / ,E 2 ) <$(0)17^7^ 1^(0)) e 4 ^ 1 "^ 4 


CT,Si,S2,<T' 


(57) 


where 

(^(0)|7^7bI^(0)) = 

Q(E-m)(m\^Ja%\m)+ 
e(-E-m)(mrnpE\m) m 

is the particle occupation number in the initial state. 

In the present problem, in particular, the initial state 
is the ground state of the pre-quench hamiltonian, i.e. 


a a E \m) 

= /?£ $(0)) =0 and therefore 



($(0)|a^a||$(0))=0, 

(59) 

and 

($(0)|^ t |$(0)) = l. 

(60) 

Substituting this initial density of occupied energy levels, 
we find 


oo —m 

AT(t)= E J dE J dE i dE 2 j dE' w* Si (E,E 1 )w 


cr,Si,S2,cr' 


It is convenient to consider also the long time av¬ 
erage of the total number of particles lim^oo N(t) = 
lim^oo t~ l f Q dt 1 N(t'). Time averaging eliminates pos¬ 
sible persistent oscillations and makes analytical treat¬ 
ment technically easier, while keeping the qualitative be¬ 
haviour at large times the same, apart from su ch oscil¬ 
latio ns. N(t) is given by the same expressions (57) and 
(61) but with e d E i~ E 2 )t re pl ace d by 


dt'e i(El ~ E2)t ' = 


ei(Ei — E2)t _ ^ 

i(Ex - E 2 )t 


(62) 


In the thermodynamic and large time limit, the above 
expressions are dominated by the poles of the overlaps 
that are close (at distance e ~ 1 /L) to the real E and 
E' axes and within the energy windows under integra¬ 
tion. Since each of the overlap factors in the integrand 
w as (E, E r ) is highly peaked around the two resonance 
poles, at E « E' and E « E' — V , and the time averaging 
factor is highly peaked at E\ ss E 2 , there are four possi¬ 
bilities to match the four energies E,Ei,E 2 and EE\ 
matches with E 2 , while E and E' independently match 
with either Ei 2 or E\ >2 — V (i.e. E' matches with either 


t (E,E 2 )w ir , sl (E , ,E 1 )wl, a (E , ,E 2 ) e ^~ E ^. (61) 


E + V or E or E — V). But integration over E and E' is 
restricted to the windows (+m, +oo) and (—oo, — to) re¬ 
spectively, so there remains only one possibility that leads 
to overlapping energy windows: E ss E\ ss E 2 ss E' + V. 
This possibility exists only for V > 2 m and constraints 
the variables E\ ss E 2 « E to be within the Klein zone 
while E' £ (m — V, — to) (Fig. pi). Therefore, using this 
approach, the expression for N(t) reduces to a single en¬ 
ergy integral over the Klein zone only (see Appendix [B|for 
details). This clearly shows that, in the thermodynamic 
and large time limit there is no steady-state particle pro¬ 
duction for V < 2m, i.e. for the values of potential which 
do not admit the existence of the Klein Zone. 

Substituting the expressions for the resonance poles of 
the overlaps, using the properties of the coefficients A, 
and Bi and performing the integration using the residue 
theorem (Appendix |b]), we finally find the simple expres¬ 
sion 


m = -t 


V—m 

J T{E)dE 


m 


(63) 
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E E u E’ 



FIG. 5: Resonance between energies in the integral (57 1 . E 
and E' match with Ui ,2 in four possible ways, indicated with 
arrows. Since they are integrated over the windows shaded in 
blue, the only allowed possibility is E « E\ « E 2 « E' + V, 
indicated with red arrows. Then £ 1,2 are restricted to values 
within the Klein zone, shaded green. The grey energy zones 
(— m, +m) and (V — m, V + m) indicate the totally reflecting 
zones. 


where the transmission coefficient T(E) in the Klein zone 
is given by ( |22] ). 

Our result shows that, if V > 2m, N(t ) increases lin¬ 
early with time at a constant rate, so the large time limit 
is described by a Non Equilibrium Steady State [27]. This 
linear increase for V > 2m is obviously a direct conse¬ 
quence of the existence in the post-quench hamiltonian of 
classically forbidden scattering in the Klein energy zone. 
Moreover, the rate of production of particles 

_ V—m 

A = lim = 2 lim [ T{E)dE, (64) 

*-► 00 t t 2 ?t J y K ’ 


Ennu in a three-dimensional system of Dirac fermions 
with a rounded step-like potential. The authors compute 
numerically the total electron population and find that, 
for a supercritical potential (V > 2m), there is a linear 
growth at constant rate for t —> 00 while for a subcritical 
potential {V < 2m) there is no such growth. Even though 
their setup is not exactly the same as ours, the authors 
of [30] find that their numerics are described by precisely 
the same formula (64) where T(E) is the corresponding 
transmission coefficient in their setup. Moreover, it can 
be easily checked by a trivial generalization of our calcu¬ 
lation that our result (64) is still valid when the initial 
value Vo of the homogeneous potential is not zero but 
between zero and V. 



E-m 

V-2m 


FIG. 6: The transmission coefficient T(E) and the reflection 
coefficient R(E) in the Klein zone (m, V — m), as a function 
of the rescaled energy (E — m)/(V — 2m), for several values 
of V/m. Note that T(E) —> 1, R(E) —> 0 for E = V/2, as 
V —> 00 . 


C. Physical interpretation 


turns out to increase when V increases, in accord with 
the behaviour of the transmission coefficient in the Klein 
zone (see fig. [6]). For V much larger than 2m, the particle 
production rate, indicated with A 00 (K), is given by 


Aoo(E) = 


V - 8m/3 
2tt 


V 2m, 


(65) 


which is linear in V. In Fig. [7] we show the particle pro¬ 
duction (continuous red line) and its asymptotic value 
(blue dashed one) as a function of the potential. The for¬ 
mer is exactly vanishing for V/m < 2 and it approaches 
its asymptotic value, as shown in the inset, with a preci¬ 
sion of order 10 -2 when V/m — 10. 

Our analytic result (64) is in agreement with recent 


numerical computations of the electron production rate 


The physical interpretation of the above results can 
be based on a semiclassical approach [MJ [2S], which has 
proved to give correct results in many cases of quenches 
including in inhomogeneous problems. In the semiclassi¬ 
cal approximation we view quantum excitations as quasi¬ 
particles with specified position and momentum that 
move ballistically with a velocity given by the group ve¬ 
locity that corresponds to their momentum v g < c = 1. A 
global quench injects an extensive amount of energy into 
the system by creating such quasiparticle excitations. 
These quasiparticle excitations are created at the time 
of quench, typically in pairs of coherent quasiparticles 
with opposite momenta (in the case of non-interacting 
systems). In the thermodynamic and large time limit 
where dephasing has typically eliminated quantum co- 
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FIG. 7: The large time rate of particle production A as a 
function of the potential V/m. Note that A ~ V as V —> oo. 
The dashed line shows the asymptotic linear increase while 
in the insect the difference between the two is reported as a 
function of the potential. 


herence effects, physical observables can be expressed as 
classical probability averages of the contributions of all 
such quasiparticles. This approximation is also valid for 
any time and system size, in the limit of low densities 
where coherence effects are unimportant. 

It is instructive to analyse, at the semiclassical level, 
what happens in a homogeneous quench of the poten¬ 
tial from zero value to V. In this case the Fermi level is 
shifted by a constant value equal to the height V of the 
potential, so that in comparison with the post-quench 
ground state, the (filled) pre-quench Fermi sea would 
correspond to a depleted or raised Fermi sea, depend¬ 
ing on the sign of V. From the structure of the overlaps 
in the homogeneous case we deduce that a pre-quench 
creation or annihilation operator can be expressed in 
terms of post-quench operators with the same momen¬ 
tum; hence the semiclassical picture that emerges here 
is that the initial state is a source of non-paired quasi¬ 
particles. We highlight that this is slightly different from 
the typical semiclassical picture where the initial state 
can be seen as source of oppositely moving pairs of co¬ 
herent quasiparticles. We would like to emphasize that, 
in the homogeneous quench case, there is no net particle 
production. What happens is that, due to the shift of the 
potential and the consequent shift of the Fermi energy, 
the pre-quench excitations may have become particle or 
antiparticle according to the post-quench Hamiltonian. 

This picture can be applied also in the case of a step¬ 
like profile of the potential, even though in this case the 
initial quasiparticle density is inhomogeneous. The sit¬ 
uation is depicted in Fig. [8j where it is shown that the 
particle density p(x, t ) is given by the sum of three con¬ 
tributions: direct particles (D) that have not passed from 
the origin before arriving at the space-time point (x,t), 


reflected particles (R) that have been reflected at the ori¬ 
gin and transmitted particles (T) that have been trans¬ 
mitted/refracted through the origin. To obtain the long 
time behaviour of observables, we can neglect the con¬ 
tribution of the quasiparticles produced in the quench 
from the region near the origin and consider just those 
produced from the regions at x —» ±oo. In these re¬ 
gions, where the potential is spatially constant, the initial 
quasiparticle production is that corresponding to an ho¬ 
mogeneous quench to the local value of the potential. Af¬ 
ter the quench the quasiparticles travel ballistically with 
group velocity v g with v g < c and, when they arrive at 
the origin, they scatter with the potential step. It is the 
scattering of the quasiparticles produced after the quench 
with the potential step that allows net particle produc¬ 
tion. If V > 2to the scattering of antiparticles incoming 
from the right produces particles with probability T(E ). 
Therefore the number of particles increases with time and 
more precisely it increases linearly, since for large times 
the number of incoming excitations reaching the origin 
is the number of all initial excitations created at large 
distances. We conclude that the semiclassical approach 
correctly describes the physics of the constant rate of 
particle production connected to Klein tunneling. 



FIG. 8: Space-time diagram representing the three possible 
origin for a quasi-particle passing at a point ( x, t) in the semi¬ 
classical picture. 


V. CONCLUSIONS 

In this work we study the Klein tunneling by perform¬ 
ing an inhomogeneous quantum quench of the potential 
in a system of relativistic one dimensional fermions. We 
determine the full energy spectrum of the inhomogeneous 
system and explain the necessity for a special definition 
of the eigenstates in the Klein zone due to the orthonor¬ 
malization requirement. We further derive the overlaps 
of these eigenstates with the homogeneous ones, which al¬ 
lows us to study the evolution of physical observables in 
the thermodynamic and large time limit, exploiting the 
resonances of pre-quench with post-quench eigenstates. 
In particular, we derive analytically the particle produc¬ 
tion rate A and show that its dependence on the ratio 
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V/m reflects the paradoxical fact that the transmission 
of relativistic fermions inside a potential step is nearly 
complete, when the step size is large. We have therefore 
demonstrated that a quantum quench of the potential 
can be used to reveal the presence of Klein tunneling, 
especially since the extensive amount of energy injected 
by the quench leads to a steady particle production. 

Our method shows that the resonance transitions are 
dominant in the thermodynamic and large time limit. 
From this point of view, it is reminiscent of the semiclas- 
sical theory of radiation and Fermi’s Golden Rule, even 
though the latter relies also on perturbation theory, while 
our method does not. This explains why our results are 
consistent with the semiclassical interpretation of quan¬ 
tum quenches, as we also showed in more detail above. 
Following this observation, we expect our result in the 
form (641 to be valid even in different spatial dimensions 
and for different profiles of the potential step. In fact, 
as we mentioned, (64) is in perfect agreement with the 
formula that was found to describe numerics in a three- 
dimensional analogue of our problem with a smoothed 

step m, 

Although the major challenge of the steepness of the 
potential is still present even in our quench setting, we 
have analytically shown that the physics of the Klein 
tunnelling can be obtained in a context very different 


from the usual one. The novelty of our approach consists 
in the fact that by quenching the value of the potential 
we no longer need to have an incident beam of particles. 
A possible experimental realization could involve highly 
controllable systems in table-top experiment with optical 
lattices and ultracold atoms. 

It would be interesting to investigate the consequences 
of quenching other parameters of the Dirac equation, as 
for example the mass or the mass phase of the relativistic 
Dirac fermions from a constant value to an inhomoge¬ 
neous kink like shape [36] . Actually vacuum polarization 
effects involving fermionic fields interacting with back¬ 
ground solitons have been shown to induce fractional 
fermion number localized on the soliton 137. Follow¬ 
ing the stimulating proposals of J3S] and [.'I9j it would 
be interesting to observe the consequences of a quantum 
quench on the charge fractionalization mechanism. 
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Appendix A: Finite volume energy eigenstates 


As typically in quantum field theory, infrared divergences occurring in observables can be resolved by setting up 
the problem on a finite size system. Momentum or energy integrals are then replaced by sums over discrete sets of 
eigenvalues, determined by the quantization conditions, i.e. the boundary conditions of the finite size system, which 
are assumed periodic. In the present problem this means that in determining the energy eigenstates, apart from the 
matching condition at the origin, we also have one at x = ±L/2, where L is the length of the system. Following the 
same procedure as before, we write the eigenstates in the general form 


,, f ) = I (E- x ) + Bu 0 (E; x ), 

1 Cuy(E; x) + Duy {E\ x), 


x < 0 
x > 0 ’ 


(Al) 


and impose the conditions 


ip(x —» 0 + , t) = ip(x —> 0 , t), 
ip(x -A —L/ 2, t) = ip(x -A +L/2, t ) 


(A2) 

(A3) 
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from which, eliminating the coefficients A, B , C and D , we find that the wavenumbers p = k,E and q = Ue-v must 
also satisfy the condition 


sin 2 [(p + q)L/4] 
sin 2 [(p- q)L/A] 


1 — K 
1 + K 


(A4) 


with k given in terms of p , q and E in ([£]). Expressed in terms of the only independent variable E , this is the equation 
for the energy eigenvalues at finite size L. Note that it is a transcendental equation and can be solved numerically, 
although several qualitative results about the solutions can be derived without explicitly solving it. We notice that due 
to the periodicity of trigonometric functions, the distance between two successive quantized momenta is of the order 
7 r IL and therefore the density of states in momentum space is uniform. Other results, like the absence of solutions 
in the energy window ( V — m, m) for V < 2 m and the degree of degeneracy of solutions, can also be concluded from 
the form of the equation. 

In evaluating expressions for physical observables, we can use a standard trick SO] in order to pass from discrete 
energy sums to integrals in the thermodynamic limit. The energy sum can be written as a sum of residues of a 
complex function with single poles at the roots of the equation (A4|. Next the sum of residues can be written as a 


contour integral, more specifically a sum of two integrals, one above and one below the real energy axis, excluding 
any other possible poles that are on or close to the real axis (at distance ~ 1/A). In the thermodynamic limit, one 
of these integrals vanishes because the integration contour can be deformed to imaginary infinity where it decays 
typically exponentially with L. If there are any singularities away from the real axis that are crossed while deforming 
the contour, their contribution vanishes also exponentially with L. The other of the two integrals gives the desired 
thermodynamic expression for the observable, along with the correct prescription of how the contour passes around 
any pole that is on or close to the real axis. 

Applying this method in our present problem, firstly justifies the regularization of infrared divergences with e ~ 1 /L 
used in the main text, since this is the order of the distance between successive energy eigenvalues. Secondly it justifies 
that any poles whose distance from the real axis is of order more than e do not matter in the thermodynamic limit 
and in fact they typically decay exponentially with L. 


Appendix B: Evaluation of energy integrals 


We will show how the asymptotics of the multiple energy integrals in (571 in the thermodynamic and large time 


limit, can be derived by deformation of the integration contours in the complex energy plane and using the analyticity 
properties of the integrated function, more specifically the resonance poles of the overlaps which we studied in Section 
1 A The following analysis would give the exact value of the multiple energy integral for all times and for large 
but finite system sizes, if these poles were the only singularities of the overlaps in the complex energy plane. This 
is not true, first because the overlaps are piecewise functions, i.e. they have a different expression in the Klein zone 
in comparison with the others, and second because the functions D ± (E,E ') 1 A{E) and B(E ) exhibit branch cuts 
and poles away from the real axis. These non-analyticities prevent the deformation of the integration contours in the 
complex energy plane. However, given that we are only interested in the limit e ~ 1/L —> 0 and t —> oo, the asymptotic 
form of the integrals is given by the contribution of the resonance poles and therefore it is sufficient to know that 
A(E) and B(E ) are simply bounded functions for all energy values and smooth at the poles. Indeed, as argued in 
Appendix |Aj the contribution of any singularity away from the real energy axis is suppressed in the thermodynamic 
limit. For the same reason the integration ranges can be extended from any energy zone to the whole real axis 
(—oo, +oo), provided this does not introduce additional resonance poles, because the asymptotics are not affected by 
this extension. It should also be noted that the piecewise form of the overlaps can be understood by expressing the 
eigenstates as functions that are analytical in a non-trivial Riemann surface with suitably chosen branch-cuts, based 
on physical requirements [25]. This suggests that it may actually be possible to deform the integration contours, so 
that our resonance pole approach that we present below can be extended to give exact results for any time, not only 
at large times. Such an exhaustive analysis is however beyond the scope of our present work and it is sufficient for us 
to consider this approach as a widely used approximation. 

The first step in evaluating the integral is to substitute the expressions for the overlaps keeping only the relevant 
energy poles, i.e. those close to the real energy axes. These poles lead to matching of the energies in pairs, so the 
next step is to identify all possible energy pairing possibilities that are consistent with the ranges of the energy zones 
under integration. We thus find E & Ei & E2 ~ E' + V , so that we can replace all energy variables by a single 
variable everywhere except at the poles, since the positions of the poles and the prescription for the integration around 
them, determines the asymptotic L and t behaviour. These substitutions are essentially the outcome of evaluating 
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the residues of the poles and can be performed directly. We finally find three types of integrals: 
' dEdE 1 dE 2 dE' 


h = 


, (27T) 4 

gi(Ei — E2)t ^ 


Fi(E)x 


h = 


i(Ei E 2 )t (E — Ei + iepp}) (E 2 — E + iep E l ) (E r + V — Ei + iep E V ) (E 2 — E' — V + itp E V ) 
dEdE-[ dE 2 dE' 

j —(sy — Fl{E)x 

gi(Bi-B 2 )t _ 1 i i i i 

i{Ei - E 2 )t (E-Ei+ iep^ 1 ) (E 2 ~E + iep^ 1 ) (Ei — E' — V + iep^y) {E' + V-E 2 + iep^y) 
T f dEdEidE 2 dE' 

h= J —(sy— Fs(E)x 

g z(Bi -E 2 )t _ 1 j i l i 

i{Ei - E 2 )t (E 1 ~E + iep^ 1 ) {E~E 2 + iep(Ei — E' — V + iep E \ v ) (E' + V-E 2 + iep^ly) 


(Bl) 


(B2) 


(B3) 


Each of the three integrals can be evaluated using the residue theorem. First we perform the integration over E and 
then over E ', closing the integration contour either above or below the real axis.Then we can perform the integration 
over Ei using the residue theorem in the form 


+oo / 

P.V. J dz f(z) = 2iir I Res f(z) + ^ Yj Res /( z ) 


real axis 


(B4) 


which is valid for any function f(z) that decays sufficiently fast for ^ —» *oo, as in the present case, where the integrand 
decays exponentially as Ei —>■ ioo. The remaining integration cannot be performed without explicit knowledge of the 
functions F t . However, provided that convergence allows it, we can postpone the integration for later and take first 
the thermodynamic limit e —> 0 followed by the large time limit T — > oo to find the asymptotic behaviour we are 
looking for. 

Following this procedure for each of the three types of integrals, we find for the above asymptotic limits 


V—m 


= 1‘ 1 t F ^ E) ’ 

(B5) 

V—m 

h = / ^F 2 (E)(-^tp-A— -^=3-, 

J 2n \e ) p E + Pe~v 

m 

(B6) 

O 

II 

(B7) 

The explicit expressions of the functions Fi(E) and F 2 (E) are 


Fi(E) = \B 2 (E)\ 2 =T(E), 

(B8) 

F 2 (E) = | Ai(E)Bl(E) + A*(E)B 2 (E)\ 2 = 0, 

(B9) 


where the coefficients Ai and Bi are given by S 
Summing up all terms, we find our final result (63) 


since the remaining energy integration is over the Klein zone. 


N(t) 


V—m 

i* j ns)iE. 


m 


(BIO) 













